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The Input: the Sequences S*^; (n) 

Recall that a Young Shape (alias Ferrers diagram) X ~ (Ai,A2,...) is an arrangement of left- 
justified rows of boxes, with Ai boxes in the first row, A2 boxes in the second row, etc. If the total 
number of boxes is n, we write A h n. 

The set of shapes A with n cells with the property that Afe+i < I is denoted by H{k, I; n). 
Recall also that a Standard Young Tableau (SYT) of shape A h n is any way of placing the integers 
from 1 to n inside the n boxes in such a way that all rows and all columns are increasing. There is 

^j , a well-known, explicit, formula for /^, the number of SYTs of a given shape A, due to Young and 

Frobenius, as well as a more elegant reformulation, called the hook-length formula, due to Frame, 
Robinson, and Thrall. 

SiJ \ In 0, Regev and Berele studied the sequences 
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and obtained explicit leading asymptotics for general (k, I) as well as general (real) z. 
Except for (fc, G {(1, 0), (2,0), (4,0), (1,1)} (when z = l)and(fc,0 S {(1,0), (2,0), (1, 1), (00, cx))} 
(for z = 2), there are no known "nice" expressions for Sffl{n) in terms of n. Of course, using the 
Young-Frobenius formula one can always express it as a multi-sum of hypergeometric terms, but, 
in general with fc + ^ — 1 sigma signs. 

\^ , The next best thing to a closed-form formula is a linear recurrence. It follows from the Funda- 

O^l ' mental Theorem of Wilf-Zeilberger theory ([3]) that for any specific {k,l) and positive integer z, 

. . the sequence A{n) ~ S)^l{n) is holonomic (alias P- recursive) , in other words it satisfies a (homo- 

C^ I geneous) linear recurrence equation with polynomial coefficients. This means that there exists a 

^Z. • positive integer L, and polynomials po{n), . . . ,p]^{n) such that 
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Furthermore there exist algorithms for computing these recurrences, but they are rather slow for 
C^ ' multi-sums. 

Since we know a priori that such a recurrence exists, we may just as well find it empirically by 
asking the computer to compute the first 60 (or whatever) terms and then guess the recurrence 
by undetermined coefficients. This is done by the command listtorec in Bruno Salvy and Paul 
Zimmermann's Maple package gfun (but we prefer to have our own version). 
While there is a "theoretical" possibility that the guessed recurrence is not the "right" one, this 
is very unlikely, and besides, using the recurrence one can extend the sequence very fast, and 
check whether the continuation matches, and if it does (say for the next 20 terms), then this is a 
semi-rigorous proof, good enough for us. We know that we have the option to find a fully rigorous 
proof by the multivariable Zeilberger algorithm [T] (that alas, will take much longer, and for larger 
{k,l) may not terminate in a reasonable amount of time) but why bother?. 
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The beauty of knowing a linear recurrence is that one can use the Poincare-Birkhoff-Trjitzinsky 
method to derive higher-order asymptotics for the sequence. This method, described in [3], is fuUy 
implemented in Doron Zeilberger's Maple package 

http: //www. math. rutgers . edu/~zeilberg/tokhniot/AsyRec 

The drawback of that method is that it does not output the constant factor, that must be deter- 
mined empirically, in general. Luckily, thanks to [2], we know that constant exactly (expressed in 
terms of tt), and this enables the computation of very refined asymptotics of the sequences Sj^ [ (n) 
for the most interesting cases z = 1 and z = 2, and for fc + ^ < 5. 

The Output: Recurrences and asymptotics for s[^j{n) and SJ^jin) for k + l <5 

All the necessary algorithms have been implemented in the self-contained Maple package HOOKER, 
kindly written by Doron Zeilberger, and available from 

http: //www. math. rutgers . edu/~zeilberg/tokhniot/HOOKER 

Detailed output files, with lots of terms in the sequences, and asymptotics to order 10, can be 
viewed in the webpage of this article: 

http: //www. math. rutgers . edu/~zeilberg/mamarim/mamcirimhtml/hooker .html 
Here we only list (for the record) the output for (fc, I) E {(2, 1), (2, 2)} and z = 1, 2, and even for 
these, we only list the first few terms, the annihilating operator and the refined asymptotics to 
order 3. 

Below A^ denotes the shift-operator in n, so, for example, the annihilating operator for the Fi- 
bonacci sequence is A^ — A — 1 and for n! is A^ — ?i — 1. 

(fc,Z) = (2,l), z = l 

1, 2, 4, 10, 26, 71, 197, 554, 1570, 4477, 12827, 36895, 106471, 308114, 893804, 2598314, 7567466, 22076405, 

64498427, 188689685 
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(fc,/) = (2,l), z = 2 



1, 2, 6, 24, 120, 695, 4403, 29540, 206244, 1483371, 10919271, 81896661, 623810421, 
4813777566, 37561178658, 295907998908, 2350767037116 

{n + 2f {l9n^ + A0n+18){n + 2fN 

{n + 3f {n + 3f{n + lf 

(45 + 148 n + 159 n2 -K 70 n^ 4- 11 n^) A^2 
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(fc,/) = (2,2), z = l 

1, 2, 4, 10, 26, 76, 232, 764, 2578, 9076, 32264, 117448, 428936, 1589680, 5897504, 22101304, 82851218, 
312935236, 1182083272, 4491680504, 17067914056, 65167445872 
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(for more refined asymptotics see the webpage). 

(fc,/) = (2,2), z = 2 

1, 2, 6, 24, 120, 720, 5040, 40320, 361116, 3540600, 37207368, 411988896, 4747167568, 
56428884512, 687793860000, 8559142303296, 108400653865572 

(2n-l)(n-l)n (231 n + 711 n^ + 164 n^ + 588 n^ - 2) N 

_5]^2 -^ — — + 16 -^ 

(n + 5)(n + 4)^ (n + 5) (n + 1) (n + 4)^ 

(1146 n^ + 545 n^ + 1003 n + 94 n^ + 296) N^ (2030 n? + 85 n" + 688 n^ + 2501 n + 996) N^ 

{n + 5) (n + 1) (n + 4)^ {n + 5) (n + 1) (n + Af 

(87 n^ + 8 «■* + 525 n + 336 71^ + 250) iV* 

(n + 5)(n+l)(n + 4)^ 

woo.^n/ -n7/2 A 33 _i 2145 _2 81723 ^oX ./^ 
l/3216"(n 1) [l + -n ^ + ^n ' + j^ n ^jn ^^ 
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